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Abstract 

A replicator equation with explicit space and global regulation is considered. This model 
provides a natural framework to follow frequencies of species that are distributed in the space. 
For this model, analogues to classical notions of the Nash equilibrium and evolutionary stable 
strategy are provided. A sufficient condition for a uniform stationary state to be a spatially 
distributed evolutionary stable state is presented and illustrated with examples. 
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1 Preliminaries and notation 

The general replicator equation comprises well-established biomathematical models that arise in 
quite distinct evolutionary contexts (see, e.g., [13 [181 [25]). I n particular, this equation appears 
in the areas of theoretical population genetics (e.g., [TTJ, [26]), prebiotic molecular evolution 
[!J U Q2], and evolutionary game theory fLH \ l22 | 127] . 

Arguably, one of the simplest replicator equations takes the form 



:av)< - f io %t) 



1 n. (1) 



Here v = v(t) = (vi(t), . . . , v n (t)) E M. n is a vector-function of n variables, A is a constant nxn 
matrix with elements aij £ M, (Av); is the i-th element of the vector Av, (Av); = Y27=i a ij v j{t), 
f loc (t) is a function, which is determined later, and dot, as usual, is used to denote differentiation 
with respect to the time variable t. It is customarily supposed that Vi(t) describes a relative 
abundance of the i-th species such that the total concentration Y17=i v i(t) ^ s kept constant and 
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often equal, without loss of generality, to 1. Therefore, the state space of (P) is the simplex 
S n = {v: Y^ 1 i=i v i(f) = 1) Vi(t) > 0, i = l,...,n}, which is invariant under (TjQ) if we set 
f loc (t) = (Av, v) = ^™ =1 (Av)ji>j, so that (•,•) denotes the usual scalar product in W 1 . Since 
(Av)j gives the net rate of growth of the i-th species in our system, the quantity (Av), is termed 
as the Malthusian fitness, and hence f loc (t) represents the average fitness of the population at 
time t. 

System ([I]) is a very well studied object, see, e.g., [8j[T8j[l7]. Very briefly, the rest points of 
(PQ) are given by the solutions to the following system 

(Av)i = (Av) 2 = . . . = (Av)„ = (Av, v) = 0, v G S n . (2) 

In general, ([2]) can have no solutions in int S n = {v : Yli=i v i = 1j v i > 0> i = 1, • • • ,n}, & unique 
solution, or infinitely many solutions. The necessary condition to have a solution to ([2]) is the 
linear dependence of the rows of A, which we denote Aj, and the vector l n = (1, 1, . . . , 1): 

n n 

MoJ-n + y^MiAj = 0, /Uo^O, y~]fH ^ 0. 

i=l i=l 

A natural way to derive the replicator equation ([1]) from the first principles is to start with 
a selection system of the form 

yi = Fi(v)yi, i = l,...,n, (3) 

where y = (yi(t), . . . , y n (t)) G R™ is a vector of absolute sizes, and i^(v) denotes the Malthusian 
fitness (the per capita birth rate) of species yi, which can depend on the structure of the total 
population at the time t. Assuming that y does not tend to zero, the change of the variables Vi = 
Vi I (J27=i Hi) l e& ds to the replicator equation ([1]), if Fj(v) = (Av)j. Therefore, it is equivalent 
to study the selection system (|3|) or the replicator equation ([TJ (see [5DJ [21] for more details) . 

As it was mentioned, the replicator equation ([Tj) naturally arises in the evolutionary game 
theory (for the origin, see [221 123] ). There exists a parallel between the concepts of the game 
theory with a payoff matrix A and the behavior of the solutions to the replicator equation ([T]) 
[17j . In particular, one of the central notions of the game theory, the Nash equilibrium v, is 
defined as such v£S„ for which 

(v, Av) < (v, Av) (4) 
for any v G S n ; and an evolutionary stable state (ESS) v G S n is defined as 

(v, Av) > (v, Av) (5) 

for any v ^ v in a neighborhood of v G S n . 

It can be shown (see |17] for the details and proofs) that if v is a Nash equilibrium of the 
game with payoff matrix A, then v is a rest point of ([T]). Moreover, if v is a rest point of ([T]) 
and Lyapunov stable, then it is a Nash equilibrium; and if v is ESS, then it is an asymptotically 
stable rest point of ([I]). 
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We remark that model ([I]) is a system of ordinary differential equations (ODE), i.e., it is 
a mean-field model. A significant attention was drawn to the replicator equation in the case 
when heterogeneous spatial structure can be included into the model formulation One of 
the suggested solution was spatially explicit models (see [H [21 [10] for the models of molecular 
evolution). In general, there are several different approaches to include spatial structure into 
the replicator equation. The solution to the problem when all the diffusion rates are equal 
is straightforward: in this case, following the ecological approach, we can just add the Laplace 
operator to the right hand sides of ([I]) (this was used, e.g., in |141I16|). However, the assumption 
of the equal diffusion rates would be too stringent in the general situation. To overcome this 
problem, Vickers et al. [H [19l [28] introduced a special form of the population regulation to 
allow for different diffusion rates. In these works a nonlinear term is used that provides a local 
regulation of the populations under question to keep the total population size constant, although 
no particular biological mechanism is known that lets individuals adapt their per capita birth 
and death rates to local circumstances [13j . A more straightforward approach, in our view, 
would be to start with a spatially explicit selection system of the form 



and apply transformation i>j = y%/(Yli=i In Vi ^ x )- Here yi = ?/i(x, t), x £ £1 C R"\ and di > 
0, i = 1, . . . , n, are diffusion coefficients. This approach is similar to the lines how the replicator 
equation (TjQ) can be obtained from ([3]). In this way, assuming impenetrable boundary of the area 
Q, we automatically obtain the condition of the global population regulation 



which was considered in some earlier works on the mathematical models of the prebiotic molec- 
ular evolution [31 [U El [29] . Here we extend this approach to the general replicator equation. 

The rest of the paper is organized as follows. In Section [2]the model formulation is presented, 
together with some additional definitions and notation. Section[3]is devoted to stability analysis 
of spatially homogeneous equilibria of the distributed replicator equation. In Section H] we 
formulate possible extensions of the notions of the Nash equilibrium and ESS for our spatially 
explicit model and present some consequences of the new definitions. In Section [5] we derive 
sufficient conditions for a distributed ESS along with some illustrative examples. 

2 Replicator equation with explicit space 

Let f2 be a bounded domain, f2 G K m , m = 1, 2, or 3, with a piecewise-smooth boundary T. 
In the following we assume, without loss of generality, that the volume of £1 is equal to 1, i.e., 
dx=l. A spatially explicit analogue to ([!]) is given by the following reaction-diffusion system 



= Fiiy)Vi + diAyi 



i = 1, . . . , n 




n 



i=l 



d t Ui = Ui [(Au)j - f sp (t)} + diAui, i 



= 1, ... ,71 



t > 0. 



(6) 
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Here U{ = Ui(x,t), x € f2, dt = , A is the Laplace operator, in Cartesian coordinates A = 
SfcLi "jj^z , A is a given constant matrix, and di > 0, i = 1, . . . , n are the diffusion coefficients. 
The initial conditions are Ui(x, 0) = (fi(x), i = 1, . . . , n, and the form of f sp (t) will be determined 
later. 

It is natural to assume that we consider closed systems (see also [29])) i-e., we have the 
boundary conditions 



o n Ui\x€r 



fin 



= 0, i>0, i = l,...,n, (7) 

xer 



where n is the normal vector to the boundary T. 

As it was discussed in Section [H the global regulation of the total species concentrations 
occurs in the system, such that 



y / «i(x, t) dx = 1 



(8) 



for any time moment t. This condition is analogous to the condition for the constant concen- 
tration in the finite-dimensional case |17j . From the boundary condition ([7|) and the integral 
invariant ([8]) the expression for the function f sp (t) follows: 

f sp (t) = E / n *( Au )* dx = f ^ Au ' u ) dx - ( 9 ) 



Suppose that for any fixed x G f2 each function Ui(x,t) is differentiable with respect to 
variable i, and belongs to the Sobolev space W^(fi) if m = 1 or Wf(^) if m = 2,3 as the 
function of x for any fixed t > 0. Here W£ , A; = 1,2 is the space of functions, which have square 
integrable derivatives with respect to i £ up to the order k. Note, that from the embedding 
theorem (e.g., [23]) we have that any function from the space (fi) is continuous, except 
possibly on the set of measure zero, and this continuity is used in some proofs below. 

Denote fit = ft x [0,oo) and consider the space of functions -B(fit) with the norm 



IvIIb = ™ax<j \\y(x,t)\\ W k + ||||(x,i;|| 1r /. 



Hereinafter we denote S n (£lt) the set of non-negative vector-functions u(x,t), Ui(x,t) £ 
B(Qt), i = 1, ... ,n, which satisfy ([8]), and we use the notation int S n (£lt) for the set of functions 
u(x, t) E S n (Qt) fo r which Uj(x, t) > 0, i = 1, . . . , n. 

We consider weak solutions to ([6]), i.e., the solutions should satisfy the integral identity 

no oo oo 

<)Ul 7] dxdt = [ ( m [(Au)j - f sp (t)} rj dxdt - di [ [ (Vui, Vt]) dxdt 



dt 

on oh oh 



for any function rj = r/(x, t) on compact support, which is differentiable on [0, oo) with respect 
to t and belongs to the Sobolev space W^fi) for any fixed t > 0. 
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Remark that generally system (|6|) is not a "system of differential equations" because its 
right-hand side contains functional Q. 

The steady state solutions to ([U]) can be found as the solutions to the following elliptic 
problem 

diAwi + Wi [(Aw), - f s P] = 0, i = l,...,n (10) 

with the boundary conditions 9 n u>j(x) = on T. Here lOj(x) G W2 (^)- 
The integral invariant ([8]) now reads 

V f w;(x)dx = 1. (11) 
i=i ^ 

The set of all non-negative vector-functions w(x) = (wi(x), . . . ,w n (x)), ifi(x) G W^ST), i = 
1, . . . ,n, that satisfy (fTTj) is denoted S n (Q). Using (QT| and (fTUj) we obtain that 

f sp = [ (Aw, w) dx, (12) 
Jn 

i.e., / sp is a constant. The rest points of JJJ) are spatially homogeneous solutions to (fTUj) ; 
the converse also holds: any spatially homogeneous solution to (jlOp is a rest point of ([1]). In 
[5j [6] we proved that for sufficiently small values of the diffusion coefficients di there exist non- 
homogeneous solutions to (fTU|) in the case (Au)j = kiUi and (Au)j = fejitj-i for arbitrary positive 
constants ki (these are autocatalytic and hypercyclic systems, respectively, for more details see 

mm)- 

Here we undertake a task to investigate the stability of spatially homogeneous solutions to 
([6]) , which can differ from the stability of the rest points of the local system ([1]) due to the explicit 
space in the model. We also consider the conditions necessary for spatially non-homogeneous 
solutions to appear. And yet another purpose is to transfer the definitions of the Nash equilibria 
and ESS for the distributed system ([6])- ([8]) and apply these concepts to study the asymptotic 
behavior of ([!]). 

3 Stability of spatially homogeneous equilibria 

We start with the standard definition that is used extensively throughout the text. 

Definition 1. A stationary solution w(x) G S n (£l) to (jlOp is Lyapunov stable if for any e > 

there exists a neighborhood 

n 

u 5 = (w(x) g s n (n) : IM X ) - Mx)\\wHti) < 5 ") 

i=l 

of w(x) that for any initial data of ([6]) from II s it follows that 

n 

Y,\\Mx,t)-Mx)\\ 2 B(n t )<z 2 (13) 

i=l 
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for any t > 0. 

If in (|13p i/ie Ze/t hand side tends to zero, then w(x) is asymptotically stable. 

In Definition [JJ Uj(x, t), i = 1, . . . ,n, are the corresponding solutions to © with the initial 
data Wi(x) G Z7" 5 , i = 1, . . . , n. 

Consider the following eigenvalue problem 

AV(x) + AV(x) = 0, xeSl, a n V| xG r = 0. (14) 

The eigenfunction system of ([T4"|) is given by f/>o( x ) = 1) {^(x)},^ and forms a complete system 
in the Sobolev space W|(0) (e.g., [21]), additionally 

(^i(x), Vj(x)) = / ipi{x)ipj(x)dx = Sij, (15) 

where 5y is the Kronecker symbol. The corresponding eigenvalues satisfy the condition 
= Ao < Ai < A2 < . . . < Aj < . . . , lim Aj = +00. 

i— >oo 

The following theorem gives a necessary condition for a spatially homogeneous solution to 
([6]) be asymptotically stable, if the corresponding rest point of ([JJ is asymptotically stable. 

Theorem 1. Let v G int S n be an asymptotically stable rest point of ([TJ. Then for this point 
to be an asymptotically stable homogeneous stationary solution to ([6]) it is necessary that 

n a 

V^>f, /3 = (Av,v), (16) 
i=i Al 

where \\ is the first non-zero eigenvalue of (|14p . 

Proof. Remark that if v is Lyapunov stable then v is a Nash equilibrium of the game with 
the payoff matrix A. Consider the solutions to ([6])-([8|) assuming that the Cauchy data are 
perturbed: 

Ui(x, 0) = (fii(x) =Vi + Wi(x), i = l,...,n, 

where <(x) G W*(Q) such that £? =1 IK( X )II^(Q) ^ S ^ 6 > °' 
Let us look for a solution to {(6j) — dHI) in the form 

00 

iti(x,t) = Vi + Wi(x,t), Wj(x, t) = cl(t) + y] cl(i)ip k (x), i = l,...,n, (17) 

k=\ 

which is always possible since the eigenfunctions t^i(x) of (|14|) form a complete system in (£2). 
c fc(^)' ^ = 1) 2, 3, . . . , i = 1, . . . , n are smooth functions of t. Note that w®(x) = Wi(x, 0), i = 
1, . . . , n. The spatially homogeneous equilibrium v is stable if c l k {t) — > 0, t — > 00 for all i and Zc. 
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Due to the fact that 

n „ n 

y2 Uj(x,t) dx = y^vj = i, 
l= i J n i=i 

we have from (|14p that 

n 

E4(t) = o. (18) 

4 = 1 

Substituting (|17p into Q and retaining in the usual way only linear terms with respect to W{ 
we obtain the following equations: 

rlr 1 (t) 

^ii = u t [(Ac )i - (A^v, c )] - fii(Av, c ) + c (t) [(Av)» - (Av, v)] , (19) 

for i = 1, . . . , n. Here Co = (cj)(t), • • • ,£$(£)), and r denotes transposition. From ([2]) it follows 
that the last term in f)19j) is zero. From Q and (|18p we also have that 



(Av, c ) = £4(t)(A*)i = /3^4(t) = 0. 



i=l i=l 

Therefore, we obtain that 

dd (t) 

— ~u~ = [(Aco)j — (A T v, Co)] , i = l,...,n. (20) 
at 

System ([20]) is linear, with the matrix Q = ||(/ij||ij=i,...,n with the elements 

Qij = aijii - (A T v)iVi, i,j = l,...,n, (21) 

which coincides with the Jacobi matrix of ([I]) at the rest point v. Due to the assumptions the 
rest point of ([1]) is asymptotically stable, therefore the trivial stationary point of (|20j) is also 
asymptotically stable. 

Multiplying equations © consequently by tp^x), k = 1,2,... and retaining only linear 
terms, after substituting (fT7|) into ©, we obtain 

dd (t) 

-^ = « i (Ac k ) < -A fc di4(t), A; = 1,2,..., i = l,...,n, (22) 

where c^ = (cjL(i), . . . , c^(i)), A^ are the eigenvalues of ()14|) . Linear system ()22[) has the matrix 
Rfc with the elements 

r ij = a ijVi ~ ><kdi5ij, 

where <5jj is the Kronecker symbol. From (|21H and the last expressions we have that 

r ij = Qij + (A T v)iVi - XkdiSij. 
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System (pO|) is stable, hence tr Q = Y17=i ?m < 0- On the other hand 

n n 

trR fc = trQ + (Av, v) - X k ^ dj < - Ai^d;, 

i=i i=i 

therefore trR^ is negative if (|16|) holds, which completes the proof. □ 

Remark 1. Theorem [1] gives only a necessary condition for the homogeneous rest point to 
be asymptotically stable. It is possible to specify the necessary and sufficient conditions for a 
particular case n = 2. In this case we have 

'an ai2 

«21 0-22 

System (TjQ) has an asymptotically stable equilibrium v (#i, £2) E int£2 if and only if 

011 < a 2 i, «22 < «12- 

Denoting c\ = an — «2i < 0, C2 = a 2 2 — 012 < 0, A = 011022 — 012021; we have 

c 2 „ Ci 

^1 = ; , V 2 



Cl + C2 Cl + C2 

System (|22|) has the matrix 

R = / an-Oi - \ k di ai 2 v\ 

k V a21 ^2 a 22«2 - Afe^y ' 

For the equilibria to be stable we need first that trR/% < 0, which is equivalent to 

J>> * (/,-- S£-). (23) 

4 = 1 x ' 

Condition ([23]) is a generalization of (fl~6j) since ciC2/(ci + C2) < 0. 
The second condition is det R& > 0, or 

d\d,2Xi — Xk(dia22V2 + e^an^i) + AV1V2 > 0. (24) 

If (^1022^2 + ^on^i) 2 > 4:did2^.viV2 then ([2l"|) is satisfied. The last inequality is equivalent to 

A > 0, ai 2 a 2 i < 0. (25) 

In the case A < the condition for the replicator equation (|6|) to have asymptotically stable 
equilibrium is that the largest root of the quadratic equation 

did 2 X 2 - \(dia 2 2V2 + 02011*1) + Avih = 
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should satisfy the condition 



A* < A 



(26) 



Therefore the necessary and sufficient conditions for the homogeneous rest point u £ intS^O) 
be asymptotically stable are the conditions ([23]) and ([25]) or the conditions ([23]) and ([26]) . 
Consider for example the simplest hypercycle equation with 



Here we have that condition (|23p is always satisfied. The condition det > can be written 
as 



If this condition does not hold then, as was shown in [6] for the general n-dimensional case, the 
inner rest point becomes unstable. 

If f2 is one-dimensional, then for the cases (Au)j = k{Ui and (Au)j = kiUi—i there exist non- 
uniform steady state solutions to ([TO]) - ([12]) (see [HE]). More precisely, a necessary condition 
for such a solution to exist is 



We note that tt 2 is the first non-zero eigenvalue of (j!4p when Q = [0,1]. Remark that this 
condition is a particular case of the general condition (j!6p . and we conjecture that the same 
situation occurs in the general case. 

The natural question is whether non-uniform stationary solutions appear in the general 
case ([6]). Here we show that, at least for some matrices A, system ([6]) possesses non- uniform 
steady state solutions. 

We rewrite (jlOp in the form 





dwi 
dx 



Pi, 



(27) 




wi [(Aw), - n , 



for any i = 1, . . . , n. The initial data are 



Pi(o) = w'M = o. 



(28) 



We assume that matrix A is such that Aw > for any w > 0. 

System (|27p is conservative, its rest points (w, p) can be found from 



Pi = 0, (Aw)j=/3, /3 = (Aw, w), i 



= 1, . . . ,n 



9 




Figure 1: 



therefore w is a rest point of ([I]). Consider the Jacobi matrix of (|27p evaluated at the rest point 
(0,w): 



where is the n x n zero matrix, I is the n-dimensional identity matrix, and J d is the Jacobi 
matrix of ([I]) at the rest point v = w when i-th row divided by d, for each row. The eigenvalues 
of J are the roots of the equation 

A 2n + det 3 l d = 0, det J l d = (di • . . . • d n )~ l det 3 l , 

where det J* is the determinant of the Jacobi matrix of (pQ) at v. if det j' < then for any 
n > 2 the last equation has pure imaginary roots. If det 3 l > then the same is possible when 
n = 2Jfe + l, k = l,2,... 

Let us introduce the following sets in the phase space: 

E = {p€R n :ft = 0, i = l,...,n}, 
U- = {weS n : (Aw),-f <0, i = l,...,n}, 
U + = {w G S n : (Aw), - fP > 0, i = 1,... ,n}, 

II = {w & S n : Wi = Vi, i = 1, . . . ,n}. 

From (I28p it follows that at the initial "time" the orbit of the system (127|) belongs to S. 
Suppose that w(0) G U~ . Then the second equation in (f27|) yields that functions pAx) increase as 
x increases, and therefore Pi(x) > 0, x > 0, i = 1, . . . , n. From the condition on A it follows that 
(Aw(0))j < (Aw(i))i, x > and therefore the values (Aw(i))j — f sp decrease as x increases, 
hence w(x) has to cross the hyperplane II where (Aw(x)); = f sp , i = 1, . . . ,n. After that, the 
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solution gets into the set U + , which implies that p%(x), i = 1, . . . , n are decreasing and w(x) are 
increasing. Therefore there exists such value x* such that Pi(x*) = w'^x*) = 0, i = 1, . . . , n (see 
Fig. QJ. 

The diffusion coefficients di characterize (n + i)— th component of the movement speed along 
the phase trajectories. An increase (decrease) in values di corresponds to the decrease (increase, 
respectively) in the movement speed of (n + i) — th component of the speed vector. Therefore, it 
is possible to find such values of di, i = 1, . . . , n such that all x* = 1. 

We remark that if we reduce values di twice this would mean that the phase orbit again 
would reach II (one cycle). Reducing di four times, we obtain the orbit that makes two cycles in 
the state space, and so on. Therefore, from the discussion above, it follows that system (fTOj) can 
have non-uniform stationary solutions which may possess an arbitrary number of oscillations 
(see a similar discussion in [6], where some examples are given). 

4 Dynamics of the distributed replicator equation 

Definition 2. We shall say that the vector function w(x) G S n (£l) is a distributed Nash equi- 
librium if 

/ (u(x,t),Aw(x))dx< / (w(x), Aw(x))dx (29) 
Jn Jn 

for any vector-function u(x, t) G S n (£lt), u(x, t) ^ w(x). 

Remark 2. Let w £ S n satisfy (|29p . Then w is a Nash equilibrium in the game with payoff A. 
Indeed, 

/ (u(x, t), Aw) cbc = (u(i), Aw), Ui{t) = / Uj(x, t)dx. 
Jn Jn 

([8]) implies that u(i) e S n for any t. Therefore we have 

(u, Aw) < (Aw, w), 
which means that w is a Nash equilibrium and therefore a rest point of ([1]). 

Theorem 2. 7/w(x) G int5 n (J7) is a Lyapunov stable stationary solution to ([6]) then w(x) is 
a distributed Nash equilibrium. 

Proof. Due to the fact that w(x) is Lyapunov stable, then, for any initial data from a neighbor- 
hood II s of w(x) in the space W£ (^), the corresponding solution satisfies (fT3j) , 

Suppose that w(x) is not a distributed Nash equilibrium. Then, using continuity of the 
scalar product, there exists an index i and constant £ > such that 

/ (Au(x,t))jdx- / (Au(x,i), u(x,i))dx > £ (30) 
Jn Jn 

for all u(x, t) G int S n (Qt) m a neighborhood of w(x) G int 5 n (i7). 
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Let u(x, t) be a solution to ([6]). Then from ([6]) we obtain 

^Et£(t) = / [(Au(x,t))i-(Au(x,t), u(x,t))]dx + di / A "; (X 'P rix, (31) 
at Jn J a Ui(x,ij 



where lnuj(i) = J*^ lnitj(x, i) dx. Using ([7]) we obtain 



r Auj(x,t) r i 



1 / <9«i(x, t) 



With the help of and (J32J), it follows from pTJ that 

^(ET^(t)) >£>0, 
at 

and hence 



2 



dx > 0. (32) 



\nui(t) > £t + ki, i = l,...,n. 
As itj(x, i) > 0, and In is a convex function, the integral Jensen's inequality implies that 



lnui(t) < lnui(t). (33) 

Therefore, 

Ui(t) >C exp{^}, * > 0, 
which is impossible because w(x) is a Lyapunov stable solution to ([6]). □ 

Definition 3. We shall say that w(x) G S'n(^) *s a distributed evolutionary stable state (DESS) 
if 

/ (w(x), Au(x,i))dx > / (u(x,t), Au(x,t)) dx (34) 
Jn Jn 

for any u(x, t) G S n (£lt) from a neighborhood o/w(x) G S n (0,), u(x, t) 7^ w(x). 

Let us introduce Definition 0J 

Definition 4. Stationary solution w(x) G S^ n (n) to (fTUj) zs stable in the sense of the mean 
integral value if for any e > £/iere exists 5 > suc/i i/iai /or £/ie initial data tpi (x) 0/ system 
([6]) ; which satisfy 

\<Pi -Wi\<5, i = l,...,n, 

where 

tfi= <pi(x)ax, Wi= u)j(x)dx, 
in 

follows that 

\Ui(t) -Wi\<£ 

for any i = 1, . . . , n and i > 0. 

i?ere = Jq Uj(x, i) dx and «j(x, t), z = 1, . . . , n are i/ie solutions to ([6]), u(x, i) G 5 n (J)t)- 
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From Definitions [T] and [J] it follows that stability in the mean integral sense is weaker than 
Lyapunov stability. For instance, consider a simple example: let g(x,t) G W^, x G [0, 1] and can 
be represented as 



k=i 

Let us suppose that co(t) — > when t — > oo. Then g(t) = g(x,t) dx — > when t — > oo, 
whereas \\g(x, t)\\ w i = X)fc=i c fcW(l + ^ 2 tt 2 ) does not necessarily tend to zero. 

Theorem 3. Let w G hit S n be a spatially homogeneous solution to (I10p (i.e., w G int S n is a 
rest point of ([T])). If w is DESS then w is an asymptotically stable solution to ([6]) in the sense 
of the mean integral value. 

Proof. Consider the set of function u(x, t) G S n (flt) belonging to a neighborhood U s of w in 
the space B(£lt). Define the functional 



We can always choose U s such that tij(x, t) > for all i because Wi > 0. Using (132j) we obtain 



oo 




n r . n 





where we suppress the dependence on t and x for simplicity. From (|3"H) we have that 



dV_ 
It 




The functional V is bounded. Indeed, applying Jensen's inequality for sums 



^Pi.lnqi < In 



i=l \i=l 



n / n 



) 



^pi = 1, qi > 



n 



we obtain 




Hence, 



v(u) = W 

l=1 Jn 



Wi lntij(x, t) (fx = Wilnui(t) < Wi In Wi. 

i=i i=i 



n n 
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bmce V(u) > for all u £ U 5 then V is a strict Laypunov functional for (|6|) and therefore 



lim hiui(t) = lnu}j, i = l, ...,n. 

t— >oo 



From the last inequality and using (J33J) we have 

lim > vbi, i = 1, . . . , n. (35) 

t—^oo 

If (|35p is strict at least for one i, then 

n n 

Elim Ui(t) > N u)j = 1 
i=l i=l 

which is impossible since 

n 

5^«i(t) = i 

for any i > 0. This proves the theorem. □ 
Remark 3. Let the conditions of Theorem [3] be met, then 

lim f sp (t) > lim f loc (t) = (w,Aw). (36) 

t— >0O t—KX) 

Indeed, in this case 

lim / (w,Au(x, t))dx= lim (w, Au) = (w, Aw). 

t— ¥00 t— >00 

Using dMD we obtain 

5 Sufficient conditions for DESS 

Application of the results obtained so far for particular distributed replicator systems reduces 
to the problem of checking conditions for DESS for spatially homogeneous stationary solution 
w G int S n 

Let us introduce the following function 

M(t)= I (u(x,i), Au(x,i))dx- / (w, Au(x, t)) dx. (37) 
Jn Jn 

If the stationary point w G int S n meets the conditions for DESS then M(t) < in a neighbor- 
hood of w in the space W^ifl) for functions u(x, t) G S n (Qt) from neighborhood of w(x) in the 
space S n (Qt), such that u(x, t) ^ w(x), x G 0, t > 0. 
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Let 

oo 

Ui {x,t) =w i + c° i (t) + J24(t)M^), i = l...,n, (38) 

s=l 

where ip s (x), s = 1,2,... are the eigenfunctions of problem (|14p . Due to (JSj) 

n 

Functions c^{t) are not equal to zero for t > simultaneously since u(i) ^ w. Additionally, 
from 

0< / u;(x,i)dx = Uiit) = Wi + c°(t) < 1, 

it follows that 

— < c°(t) < 1 — Wi, i = 1, . . . , n, 
i.e., these functions are bounded, which implies that there exists 5 > such that 

n 

\c°(t)\ 2 = Y,\c°(t)\ 2 >6 2 >0. (40) 
i=i 

Let us fix e > and consider an e-neighborhood of point w in Wj (^)- From the embedding 
theorems 

I|U(X, t) - W\\1 2 (Q) < ^ll U ( X > *) - ™\\wk(Q) < K ^ 2 ^ 

where K, K\ are positive constants. This means that 

oo 

||u(x, t) - w||| 2(fi) = (c°(i)) 2 + ]>>^)) 2 < K x e\ (41) 

s=l 

Theorem 4. The stationary point w G int is a DESS if 

(Ac°(t),c (t))<-7 2 |c°(t)| 2 , 7 >0 (42) 

/or any c°(i) satisfying (j39l) 

Proof. After inserting (|38|) into (|37|) we obtain 

oo 

M(t) = <Ac°(t),c°(t)> + E(Ac s (t),c^(t)) - (c°, Aw). 

Since w £ 5 n is a stationary state, then 

(Aw)i = ... = (Aw) n = j3, 
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therefore, using (j39|) 



n n 



(c°, Aw) = ^(Aw) !C ? = /3^>°(i) = °- 



On the other hand 



\(Ac s (t),c s (t))\ <K 2 \c s (t)\ 2 , K 2 >0, 



and, using (|42p . we obtain 



oo 



M(i) < - 7 2 |c°(t)| 2 + iT 2 |c s (i)| 2 . 



s=l 



Using ()40p and (|4ip we arrive at the estimate 



M(t) < (-7 2 5 2 + i^iK 2 e 2 ). 



If e < 7<Vv / ^i~^2, then M(t) < 0. 



□ 



Remark 4. In the case of replicator equations with symmetric A condition (|42p means that 
each orbit converges to the stationary point w € int S n [TT] . 

Example 1. Consider system ([6]) with the matrix 



The diffusion coefficients are d\ = 0.03, d 2 = 0.02. In this case (fl~6j) is not satisfied and there are 
spatially heterogeneous stationary solutions to the system ([6]). The dynamics of the solutions is 
shown in Fig. [2]^a), (b). Spatially non- uniform solutions are stable in this case. 

In Fig. [3] all possible stable spatially non-uniform solutions are presented, each of which has 
its own basin of attraction. 

Example 2. 




a < c, d < b. 



We have 
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Figure 2: Solution to ©, see Example 1 and text for details, (a), (6) Time dependent solutions; 
(c) Time evolution of f sp (t); (d) Ui(t) = L Ui(x,t) dx, i = 1,2 are shown 
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In this case, 



and flUD holds if \i < 0. 



(Ac (t),c°(t)) = ^\c°\\ 



Example 3 (Hypercycle equation). Consider a hypercycle system with three members and the 
matrix 

/ 1 
A= 10 

\ o i o 

Prom ([39]) c? = -(c§ + c|), therefore 

(Ac°(t), c°(i)) = -(c° + cf )c° - c°(c° + 4) + c°c° 

= -[(^) 2 + c^ + (c°) 2 ]<-i[(c°) 2 + (c°) 2 ], 



which yields (pE2j) . 
Example 4. 



a 


6 


c 


c 


a 




6 


c 


a 



u>i(x) 



(6) 



0.2 0.4 0.6 O.f 

X 



Wi{x) 



(d) 



Figure 3: Asymptotically stable stationary non-uniform solutions to the problem as in Exam- 
ple 1. See text for details. The numerical scheme is presented in [5] 
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Matrix A is circulant, and its eigenvalues are 

Ai = a + b + c, A 2 , 3 = (a - (6 + c)/2) ± iV3(b + c)/2. 

Assume a < (6 + c)/2, 6 > a > c > 0, be > a 2 , then Re A2,3 < 0. Ai has the eigenvector (1, 1, 1). 
This vector is orthogonal to any c° that satisfies (|39j) . hence A is negatively determined on this 
set. 

Example 5. 



A = 



/ 





a i 


«2 


0,3 


\ 









Ol 


a 2 






«2 


0,3 





«1 




V 


a i 


a 2 


03 





/ 



Let a\ + 03 > 02 . The eigenvalues of A are given by 

Ai = Oi + 02 + 02 

A2 = — oi —02 — 03 
A3,4 = -a,2 ± i(ai + a 3 ) 

In the same vein as in Example [U Ai corresponds to the eigenvector u = (1,1,1,1), which is 
orthogonal to any c° satisfying (|39f) . The other eigenvalues have negative real parts, therefore 
A is negatively determined. 

Consider the following matrix 



A 



f 


0.5 


1 


0.8 


\ 


0.8 





0.5 


1 




1 


0.8 





0.5 




V 0.5 


1 


0.8 





/ 



and the diffusion coefficients d = (0.03,0.02,0.03,0.02). Again we have spatially non-uniform 
solutions which are stable in the mean integral sense (see Fig. 0]). In Fig. [5] possible asymptoti- 
cally stable non-uniform solutions are shown. The details of the numerical scheme are given in 
®- 

Here the conditions for the spatially homogeneous solutions to be DESS are fulfilled, there- 
fore, this solution is stable in the sense of the mean integral value, which can be seen from the 
pictures. 

Acknowledgments: The research is supported in part by the Russian Foundation for Basic 
Research grant # 10-01-00374. ASN is supported by the grant to young researches from Moscow 
State University of Railway Engineering. 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

x x 



Figure 5: Asymptotically stable stationary non- uniform solutions to the problem of Example 5. 
See text for details. The solutions are such that w\{x) = w^{x) (dashed curve), and W2(x) = 
w±{x) (solid curve), f sp = 1.54 
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